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ABSTRACT 

Motivated by Harish- Chandra theory, we construct, starting from 
a simple CDD-pole ^-matrix, a hierarchy of new S'-matrices in- 
volving ever "higher" (in the sense of Barnes) gamma functions. 
These new S'-matrices correspond to scattering of excitations in 
ever more complex integrable models. From each of these models, 
new ones are obtained either by "g-deformation" , or by consider- 
ing the Selberg-type Euler products of which they represent the 
"infinite place". A hierarchic array of integrable models is thus 
obtained. A remarkable diagonal link in this array is established. 
Though many entries in this array correspond to familiar inte- 
grable models, the array also leads to new models. In setting up 
this array we were led to new results on the g-gamma function 
and on the g-deformed Bloch-Wigner function. 
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1. INTRODUCTION 

The factorized S'-matrices, which describe the scattering of excitations, 
are known for many (1 + l)-dimensional integrable models (e.g. XXZ, XYZ 
models and their generalizations [1], continuum sine-Gordon theory [2], etc.). 
The Jost functions for the corresponding input 4-point amplitudes have been 
identified [3], [4] with the explicitly calculable Harish- Chandra c-functions of 
certain, in general quantum symmetric spaces. 

Here we wish to draw attention to unexpected connections between vari- 
ous integrable models, which arise when certain coupling parameters assume 
special rational values. These integrable models then naturally arrange them- 
selves into an "array of models" . By moving from one site of this array to 
one of its nearest neighbors, the just-mentioned connections between known 
models come into full play. Sooner or later though, such moves to nearest 
neighbors in this "array of models" leads us outside the realm of the known 
and suggest the existence of new integrable models. The hierarchic array is 
presented in figures 1 and 2, and it leaves, as we shall see, a number of open 
questions. The very existence of this array is our main result. It implies 
hierarchical relations between different integrable models, which look both 
promising and instructive. 

2. SINE-GORDON THEORY AT CERTAIN RATIONAL 
VALUES OF THE COUPLING CONSTANT 

To begin with, consider soliton-soliton scattering in the sine-Gordon model. 
Up to a finite number of, for us irrelevant factors, the Jost function for this 
process is given by the infinite product [2] 
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where u is the relative rapidity, 
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and (3 is the sine-Gordon coupling constant. Now fix 
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Then 

A(u) = K ffyr , (3a) 

with a constant and 

r 2 (-) 

where r 2 is the "higher" gamma-function [5] related to Barnes' G-function 
[6] by 

T 2 (x) = (27c)-^ 2 e~ x2/2 G(x-l) . (4) 

The equations (3) are remarkable on two counts: 

i) the infinite product of eq. (la) has disappeared (it is "hidden" in T 2 ), 
and 

ii) the function c 2 which has appeared, is strongly reminiscent (see eq. 
(3b)) of the Harish- Chandra c-function [7] of the real hyperbolic plane 
(SL(2, R)/SO(2)) 

r(-) 

c(x) = nmk+k) (5) 

The difference between c(x) and C2{x) is visible from eqs. (5) and (3b): the 
gamma function T(x) is replaced by the "higher" gamma function T 2 (x). A 
similar r — > T 2 transition occurs not only for (3 2 /8tt = 2/3, Eq. (2a), but 
more generally, also for 

£- = 1--, ^ = m-l m = 4,5,6... (26) 
87T m 8n 

These are the well-known rational values of j3 2 /8n at which the spectrum 
truncates, as noted by Leclair [2]. 

The sine-Gordon soliton-soliton S'-matrix is essentially [4] the q — > 1 
limit of the S'-matrix for the scattering of two (dressed) excitations in 
the XYZ model. In terms of the parameter / introduced in [4], the regime 
(2) corresponds to the rational values 

1 = —^— m = 3,4,5,... (6) 
m — 1 
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These values of / are special for the transformation [4] 

/ -> V = - n ^n' = \ = m-l (7) 

which takes one from the magnetic model M. n (corresponding to the Z n - 
Baxter model) with parameters q, I to the M. n < model with parameters q, I', 
while leaving the S'-matrix unchanged. It is for and only for the values (6) 
of / that, starting from physical values of I and n (i.e. integer n > 2), the 
transformation (7) yields a new physical pair of values I' and n' (n' also 
integer). As will be argued in Section 5 the symmetric space underlying 
this argument is the double coseting of a Kac-Moody symmetry and the 
transformation (7) corresponds to the familiar rank-level duality [8]. The 
XYZ-like M n models yield, for I given by (6), an S'-matrix similar to (3), 
but involving a g-deformed T 2 function. The situation is similar to "S- 
wave" scattering on ordinary rank one symmetric spaces [9]. There the Jost 
function is given [10] by the ratio of two ordinary gamma functions (eq. (5) 
above), which when g-deformed is still meaningful: it corresponds [4] to the 
scattering of two first-level excitations in the n — > oo limit of the M. n model 
(the i?-matrix in this limit has been obtained recently by Shibukawa and 
Ueno [11]), or in an alternative interpretation [3], [4] to the scattering of two 
dressed excitations in the spin 1/2 XXZ chain. There is much more to this 
analogy. To understand this, we first present some old, as well as some new, 
results on higher gamma functions. 

3. ORDINARY AND HIGHER GAMMA FUNCTIONS AND 
THEIR q-DEFORMATIONS. THE q-BLOCH-WIGNER FUNC- 
TION 

The Weierstrass product for the ordinary gamma function is [12]: 

oo 

F(x) = -y{x) 1] (x + a^Waix) (8a) 

o=l 

with prefactor 

7 (x) = -e~ c \ (86) 
x 

C = Euler-Mascheroni constant, and Weierstrass factors 

w a (x) = ae x/a (8c) 
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The gamma function has half the poles of the function n/ shi7ra;, to which it 
is related by 

r(x)r(i - x) = -A- = - TT (x + a)-^ (9) 

srnra i 

(Z* = set of nonvanishing integers). T(x) obeys the functional equation 

T(x + l)=xT(x) (10) 

as well as its higher "multiplication" generalization due to Gauss. The rela- 
tions (8)-(ll) allow for two types of generalization relevant for us. 

A) The higher gamma functions 

With the notations 

r(x) = G 1 (x), x = G (x) (11) 
equation (10) can be written as 

G 1 (x + l) = G (x)G l (x) (12) 
which invites the recursive generalization 

G n+1 (x + 1) = G n (x)G 1 (x) (13) 

These equations, along with some other requirements (viz. Go(x) = x , 
G n (l) = 1, existence of all derivatives of G n and the condition that the 
(n + l) st derivative of log G n (x) should be non-negative for real x > 1 [13]) 
defines a unique sequence of meromorphic functions G n (x) n = 0, 1,2, . . . 
G2(x) is Barnes' G-function [6]. Let us call n the height of the higher G- 
-function G n . Nowadays it is customary to work not directly with the G n 
functions themselves, but with the higher gamma functions T n which differ 
from the G n by some prefactors [5]. For instance T 2 (x) is related to G 2 = G 
via eq. (4). We shall also say that the higher gamma function r„ has height 
n. 

Like the ordinary gamma functions, the functions T n also admit Weier- 
strass product representations of the form 



T n (x)= ln (x)l[(x + ar- 1 w^(x) (14) 

a=l 



5 



with 7„, given in ref. [5]. Notice that (r n (x)) _1 has poles at the negative 
integers. The pole at x — —a being of order a n ~ 1 . We can rewrite eq. (15) 
as an n-fold infinite product 



oo oo 



r„(x) 



ln{x) I] II II II (a: + a 1 + a 2 + ... + aO<L r >) (15) 

ai=0 a2=0 a n _i=0 a n =l 



The factors 7„, w^\ an can be found from the references [5], [6]. Yet another 
form of the infinite products involves intermediate lower gamma functions 
instead of x + a\ + . . . + a n = T (x + a\ + . . . a n ). For instance, we can write 



T 2 {x)=^{x)\{[T 1 {x + a)}- 1 w^\x) 



(16) 



a=l 



More generally T n is an (n — fc)-fold infinite product of functions or their 
inverses, of course with appropriate prefactors and Weierstrass factors. 

Just as the ordinary gamma function is "half" of a trigonometric function, 
eq. (9), so the higher gamma function T 2 (x), T 3 (x) . . . are also related to 
"higher" trigonometric functions A 2 (x),A 3 (x), . . . Eq. (9) generalizes to 



A n (x) = T n (x) [T n (-x)] 



cxp 



—71 



t n ~ X COt 7Tt dt 



n = 2,3, . 



(17) 

Notice the alternation of product and ratio of r n 's in the definition of A n . 
For n = l,Ai(x) is defined [5] as 7rexp[— n Jq cotnt dt] = 7r(sin nx)^ 1 which 
according to Eq. (9) is equal to r(x)T(l — x). The higher trigonometric 
function A n (x) is meromorphic of order n and obeys [5] 



Tin) 



l0g An{X) = (27Ti)"- 1 + 

where ((n) is the Riemann zeta function and 



^ (n-k)\ 



2n\ 



(18) 



(19) 



is closely related, as noted by Rovinskii [5], to the Bloch-Wigner-Ramakrish- 
nan function 



D n (y) = Re 



\k=l 



(n-k)\ 



Life(y) + 



(iogM) ? 

2n\ 



(20) 
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The function Z^y), the original Bloch-Wigner function, appears in the cal- 
culation of volumes in 3-d hyperbolic (Bolyai-Lobachewsky) geometry [15]. 
The Euler polylogarithm Li k (x), which appears in eqs. (19) and (20) is de- 
fined as 

oo n 

LU(-) = E ~ k (21) 

71=1 " 

It will be important for us that the higher gamma functions play the same 
role in the functional equation of Selberg zeta functions [16], that the ordi- 
nary gamma functions plays in the functional equation of the Riemann zeta 
function [17]. Selberg zeta functions are associated to double coset spaces 
M = F\G/K, rank G/K = 1, where K is the maximal compact subgroup 
of G and T is a suitable discrete subgroup of G. To these data one asso- 
ciates, in the simplest case of the trivial representation of T, the Selberg zeta 
function Zm(s), defined as usual, by an Euler product. For instance, for 
G = SO(2, 1), K = SO(2), this Euler product is [16], [18] 

OO 

Z M {z) = I] II (l - e-^ s+fe )) . (22) 
P ePk=o 

Here P is the set of all primitive closed geodesies on M (a closed geodesic is 
primitive if it is not the multiple traversal of another closed geodesic) and r(p) 
is the length of the primitive closed geodesic p. Unlike the prime ideals over 
which the (simple) Euler product is taken in the cases of the Riemann and 
Dedekind zeta functions, for the Selberg zeta function the product runs over 
"primitive elements" (geodesies) and is a double product. Like the Riemann 
zeta function, the Selberg zeta function also obeys a functional equation. For 
instance, for G = SO(2n, 1), K = SO(2n) [19], if we define [13], [18], [20] 

Z M (s) = Z 00 (s)Z M (s) 

zm = [r 2n ( s )r 2n ( s + i)](-^voim (23) 

(Zm as in eq. (22)), then the functional equation takes the simple form 

Z M (2n-l-s) = Z M (s). (24) 
This parallels the Riemann case, where 

C(s) = Coo(s)C(s) Us) = ^ 2 V( S -) (25) 
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and the functional equation reads 

C(l -*) = &)■ (26) 

The second relevant generalization of this "gamma hierarchy" concerns 
g-deformations. 

B) q-deformed gamma, higher gamma and Bloch-Wigner functions 



The g-deformations of G (x) = T (x) = x and of G 1 (x) = T 1 (x) = T(x) 
are classical [21] 

r , q (x) = [x] q = (27a) 

and 

oo 

r 1>s (x) = (1 - qf~ x II [1 + a] q [x + a]" 1 . (276) 

a=0 

Notice the simplicity of (27b), as compared to eq. (8). The g-deformation of 
eq. (9) we find to be 

r M 0r)r li? (l -x)= i{\ - q)(q; q)lq^ \ (28a) 

where 

q = e l27TT 

oo 

(?;?)oo = IK 1 -? ) ( 28& ) 

o=l 

oo 

7?i(u|r) = 2g 1/8 shiM - 2q n cos2u + q 2n )(l - q n ). 

n=l 

In the Jacobi theta function $ 1 we have used q related to r as in (28b) rather 
than the often encountered q = e tnr . Askey [22] has already pointed out that 
Ti, q (x) Ti >q (l — x) should involve Jacobi theta functions. Our equation (28) 
makes this fact explicit. The Jacobi $i-function is thus the g-deformation of 
the trigonometric sine function. 
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Higher gamma and trigonometric functions can be similarly g-deformed. 
For instance, we define the q- deformed higher trigonometric function as (see 
eq. (17)) 



a 2>9 (x) = n 



k=i 



T hq (x + k) 



l-q) 



-2x 



(29) 



The logarithm of A 2 will then correspond to a g-deformed T> 2 function (Eq. (19)), 
itself closely related to a g-deformed Bloch-Wigner function D 2q . We shall 
refer to \ogh 2 ^ q as a g-deformed Bloch-Wigner function, though, strictly 
speaking, it is V 2 and not D 2 whose g-deformation this provides. A similar 
g-deformed Bloch-Wigner function was considered in a different approach by 
Bloch and Zagier [23]. For the g-analogue of the n = 2 case of eqs. (18), (19) 
we then readily find 



io g A 2 »=-^ r^jiog 

1 — q Jo z { 



iq 



-1/83-1/2 fi ogz 
-vi I — — It 



2i 



(30) 



with r, (g; g)^ and $1 as in (28b) and J d q z standing for Jackson q- integration 
[21] . The g-integral of z~ x log $1 is thus the g-deformation of the dilogarithmic 
expression of eq. (19). Similar considerations apply to the g-deformation of 
the higher trigonometric functions and of the higher polylogarithms. 

4. THE ARRAY OF GAMMA AND ZETA FUNCTIONS AND 
THEIR SYMMETRIC SPACES 

With the principals all in place, we can now arrange them in an array from 
which their relations become clear. In this section we present the array as a 
mathematical construct. In the next section we shall show that to each "site" 
in this array corresponds a physical scattering problem with "geometric" Jost 
function. The array is presented in Fig. 1. 

The pivotal column of this array is the third column. Each of the entries 
in this column is a gamma function, the height of which increases in unit steps 
as we move down the column. The entry at height one (second row) is r 1? 
the building block a la Harish-Chandra-Bhamamurthy-Gindikin-Karpelevic 
of the c-function of an ordinary symmetric space (see eq. (5) above and 
ref. [7]). The function T\ has this property because tt~ s ^ 2 Ti(s/2) is the 
gamma factor (i.e. the local factor Coo(s) at the infinite place [10], [24]) in 
the function ((s), (eq. (25)), in terms of which the functional equation for 
the Riemann zeta function takes the simple form (26). 
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Descending one step down the third column lands us at the function T2, 
itself, as we have seen, the gamma factor for the Selberg zeta function Zm{s) 
of a 2-dimensional Riemann surface M 2 . Further down at height 2n we 
encounter the gamma factors of Selberg zeta functions of the 2n dimensional 
domains M 2n = T\SO(2n, l)/SO(2n). Each of these zeta functions is an 
Euler product, and the corresponding Euler factor at primitive element p (at 
height 1 this primitive element is a prime number, at height 2 a primitive 
geodesic, ...) is displayed in the second column. In the first column we display 
the Euler product, gamma factor included, i.e. what we called C( s )> Zm{s) 
in section 2. Finally, in the fourth column we present the g-deformation of 
the gamma function which appears in the third column at the same height 
(i.e. in the same row). 

The new feature in all this is, that the Euler factor (in the second column) 
at height h, coincides, up to a constant factor, with the inverse of the q- 
-deformed gamma function (in the fourth column) at height h — 1, provided 
the deformation parameter q at height h — 1 is set equal to the inverse of 
the "norm" of the primitive element at height h. (By norm of a primitive 
element we mean the prime number p itself at height h — 1, exp(r(p)) for 
a primitive closed geodesic p of length r(p) at height h = 2, etc..) This 
is shown in fig. 1 by the diagonal arrows. These diagonal arrows express a 
generalization of the p-adics-quantum-group connection [25] , [26] , [3] , [4] . 

In the array we only showed explicitly heights 0, 1, 2. At larger heights the 
same pattern is followed, with the remark that for odd dimensional spaces, 
the vanishing Euler characteristic leads to a simplified functional equation 
[16], [19], [20]. 

We can now translate all this into statements about integrable models. 

5. THE ARRAY OF INTEGRABLE MODELS 

Starting again at height 1 in the pivotal third column, we encounter the 
building block of the Jost function for ".S-wave" scattering in the hyperbolic 
plane SL(2, H)/SO(2), or equivalently for the scattering of dressed excita- 
tions in the Heisenberg XXX chain. Next to it (in the second column) is its 
adelic partner ( p (x), the building block of the Jost function 

* w = Jri) (31) 
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for "^-wave" scattering on the j9-adic hyperbolic plane SL(2, Q P )/SL(2, Z p ) 
[10], i.e. the Bethe lattice with incidence number p + 1. The Euler product 
builds the "S'-wave" scattering on the adelic hyperbolic plane SL(2, A)/ 
SL(2,Z p ) [10], or the closely related scattering on the fundamental domain 
of SL(2,Z) on the real hyperbolic plane [27]. We see thus that height 1 
(second row) in the array corresponds to SL(2), at least as far as the first 
three columns are concerned. The g-deformation in the fourth column is the 
Macdonald-Harish-Chandra [25] c-function for root system A\, again SL(2). 
This, as we saw in our earlier work [4], is modeled in the scattering of two 
(first level) excitations in the n — > oo limit of the M. n model, as well as 
[3] in the S'-wave scattering on quantum hyperboloids, or equivalently, in 
the dressed excitation scattering in the spin 1/2 XXZ model. So the full 
height 1 row in the array can be interpreted in terms of ordinary, or as in the 
fourth column, quantum symmetric spaces and modeled on specific physical 
systems. 

Now let us move down to height 2. In the third column the entry is now 
a T 2 function. As we saw in eqs. (3) this T 2 builds - much in the way that 
Ti did at height 1 - the Jost function for soliton-soliton scattering in the 
sine-Gordon model at the values (2) of the coupling parameter. What is the 
symmetric space here? The r 2 function is an infinite product of Y\ factors. 
This suggests a Gindikin-Karpelevic type product [7] over infinitely many 
roots, therefore a Kac-Moody symmetry A\. The level k is related to the 
parameter /. Specifically k = When eqs. (2a) or (6) are imposed this 
yields k = m — 1 (remember, m — \ is a positive integer). The transformation 
(7) then interchanges l~ x = k = level with n = rank +1. In other words, it 
becomes rank-level duality [8] as already mentioned in section 2. The level k 
is also related to an SU(2)q with deformation parameter q, a root of unity 
~k+2 _ ^ [28]. We are dealing with the c-function obtained from the zonal 
spherical functions of this symmetric space. Though the symmetric space 
itself is infinite-dimensional, it is truncated to a single "radial" dimension by 
the double coseting involved in going to zonal spherical functions. 

From this connection with 1-loop Kac-Moody algebras it also follows that 
the height h (of a gamma function) introduced above can be interpreted as 
a "loop number" . More specifically, the usual symmetric space of height 1 
corresponds to zero loop number, the height 2 Kac-Moody symmetric space 
just mentioned to loop number 1 (1-loop algebras). Height 3 is then to be 
associated with "2- loop-algebras" , etc... Height zero corresponds to loop 



11 



number -1. 

The T 2 function, as we saw, can be g-deformed and the corresponding 
r 2i9 function (fourth column) "builds" the Jost function for the scattering of 
dressed excitations in the XYZ model, the true n = 2 case, for the special 
values (6) of the parameter I (which as in ref. [4] is given by / = with 7 
the anisotopy parameter, and r the modular parameter)l. 

Moving now to the second column at height two, we find the Euler part- 
ners of T 2 in the Selberg zeta function Zm for which T 2 provides the gamma 
factor Z^. Again, the diagonal connection in the array shows that this yields 
the dressed excitation scattering in the XXZ model (or more generally in the 
Aioo model) as in the fourth column at height 1. Finally, at the "adelic-like" 
position (first column) at height 2 (third row) we encounter an S-matrix built 
out of Selberg zeta-functions Zm the same way that the adelic S-matrix at 
height 1 was built out of ( functions. The S'-matrix build out of Zm factors 
is new and its physical realization is not known to us. 

We can obviously continue this way, down the array and encounter higher 
(n > 2) gamma functions and their g-deformations, Euler partners and Zm 
functions. They correspond to rank one ordinary or quantum symmetric 
spaces of higher dimensionality. There is the possibility that the correspond- 
ing physical models are also higher dimensional. This would be very inter- 
esting. 

At this point, however, we prefer to take two steps back and move up to 
height zero (first row). Here r (x) = x is the building block of a Blaschke- 
CDD factor 

CBiaschkc-CDD(^) = F (ik - a) — ik - a (32) 
whence the S'-matrix 

T {ik - a) 

•JBlaschkc-CDD — T, r ■ (o6) 

r (—ik — a) 

By g-deforming it we get the type of factor we encountered at height 1 in 
the second column ("S-wave" scattering on Bethe lattice). This is again the 
diagonal connection in the array at work. Just as in Section 3, the Euler 
partners and adelic counterparts of the Blaschke factor, if at all meaningful, 
are not known to us. We are now in measure to produce the array of integrable 
models. It is presented in Fig. 2. 

This array is the main result of this section. As was emphasized, the array 
pivots around its third column. Moving down this column means enlarging 
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the symmetric spaces underlying these models. Moving to the right (to the 
fourth column) produces spin chains by q- deformation. Moving left (to the 
second column) yields "Euler partners" of the pivotal model of the same 
height, and finally (first column) "adelic-like" models. Already at height 2 
these are new models. Similarly the height 3, 4, ... models are yet to be 
physically identified. But they are there! 

6. CONCLUSIONS 

We have found here a hierarchical array of gamma functions and, in effect, 
of symmetric spaces. A symmetric space can be g-deformed (move from the 
third to the fourth column) or extended while maintaining rank 1 (move 
down the ordinary third column, or down the quantum fourth column). One 
can consider its Euler partners (move from the third to the second column) 
and its "adelic-like" counterpart (move from the second to the first column). 
To each site in the array we found a corresponding physical model at special 
rational values of some parameter. Pushing things to the extreme, we were 
led to ^-matrices for which the underlying physical model is as yet unknown: 
those involving Selberg zeta functions and those at height h > 2. 

The mathematical novelties we encountered were: the g-deformed Bloch- 
Wigner function(eqs. (29), (30)), the new g-deformed trigonometric function 
eq. (28), and most importantly the array of gamma functions and symmetric 
spaces of section 2. 
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Center at the University of Chicago for its hospitality and support. 
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